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We study small conducting particles and thin films in an oscillating longitudinal electric field. 
We find the charge, current, and field distribution in the particle, the polarizability and the electric 
dipole absorption. We account for Thomas-Fermi screening by adding a Fick's diffusion term to 
Ohm's law. Alternatively, we describe a particle as a dielectric body with a non-local dielectric 
constant which is derived in a microscopic linear-response theory. We show that both approaches 
are equivalent. 
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I. INTRODUCTION 



The classical theory of polarizability and absorption of small metallic particles in a time-varying electromagnetic 
. ■ field is a well studied and documented subject |l], JE[. The corresponding quantum theory has also been studied, 
starting with the paper of Gor'kov and Eliashberg fa] (GE). The quantum treatment can, in principle, address the 
features of (i) nonlocality of the current-field relationship, (ii) discreteness and statistics of the energy levels and (iii) 
self-consistent screening of the transition matrix elements. GE assumed that the spectral properties and the matrix 
elements were statistically independent and calculated the latter 'semiclassically' for both pure (£ > a) and 'dirty' 
(I < a) small particles. They did not take into account the screening of the external field. This was treated later by 
Lushnikov, Simonov, and Maksimenko [0 along the same lines as GE. 

ii - For dirty systems in particular, the problem has been revisited recently by Blanter and Mirlin l| (BM). They 

applied the formalism of the supersymmetric sigma model developed by Efetov H for describing the level statistics 
and the diffusive nature of the matrix elements of disordered wavefunctions (see also jjj ) . Their work confirmed the 
earlier result of GE for unscreened response, thereby justifying GE ansatz of statistical independence. In a later paper 
|3| BM incorporated screening (at zero frequency only) through an approximate formulation of the RPA in which the 

_f~. ' irreducible polarization part was restricted to first order in the diffusion term. 

q^ , Mehlig and Wilkinson [gj have given a phenomenological analysis of the electric response of small dirty systems. 

f^i ■ They considered the response to lowest order in frequency postulating in advance features of the in-phase and out-of- 
phase parts of the current relative to the electric field, namely, that they are ohmic and diffusive, respectively. (We 

^^ i will see later on that this separation of the current by phase is neither necessary nor correct in detail.) In a second 
paper [p I they restated their earlier results, starting with the nonlocal conductivity a (r, r';u>) first derived by Serota 
et. al. [||. As in the diagrammatic method of Ref. ||, the self-consistency equations for the charge/current/field 
distribution inside the particle are cast in integral equation form, which are cumbersome to apply in finite geometries 

— JL , and have not been solved. 

In this paper we show that a considerable simplification in the application of the scmiclassical theory in the 
above papers to specific geometries is achieved by modifying the classical Maxwell method for fields in dielectric 
and conducting bodies to include nonlocality in the electric current-electric field constitutive relation. We add to 
Ohm's law in the standard Rayleigh-Drude (RD) theory a Fick's law diffusion term, obtaining what we refer to 
as the Einstein transport equation. We then follow Landau and Lifshitz's observation (see G|) that an electrically 
isolated linear-response system can always be treated as a dielectric body - in this case one with a nonlocal dielectric 
response. However, it will be shown that even this response function can be avoided. Instead, after formulating 
appropriate boundary conditions, we proceed to solve Maxwell's differential equations for the charges and fields in an 
almost textbook-like fashion. Thereby self-consistency is achieved at all frequencies without the cumbersome integral 
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equations of diagrammatic schemes. The classical treatment of dielectric objects in external fields is, after all, a self- 
consistent solution method for charges and their fields. The great simplicity of this approach relative to the earlier 
work mentioned should not distract attention from the fact that it yields without analytic complication the first exact 
charge/current distributions for slab and spherical geometries and that the method is obviously applicable to more 
complex geometries. 

In Section II we discuss the Maxwell theory with the modified constitutive equation. The method of solution is 
presented in Section III, where it is illustrated for a slab and sphere. Section IV has the discussion of the equivalence 
to the earlier nonlocal integral equation first derived by microscopic theory. In Section V we compare the electric and 
magnetic dipole absorption. 

II. THE "DIELECTRIC" EQUATIONS AND BOUNDARY CONDITIONS. 

In this Section we present the modified RD macroscopic equations for a dirty metal or semiconductor which 
includes charge diffusion. The particle is assumed to be in a uniform quasistatic electric field E = zEoe~ tuJt such as 
that between parallel condenser plates. This approximates the electric response to low frequency light where b « 1; 
a is the characteristic small dimension of the particle. (The magnetic response is discussed in Sec. V.) The Maxwell 
equation for the longitudinal field and the continuity condition are 

A-Kp{r;uj) = V-E(r;w) (1) 

iwp(r;uj) = V ■ j(r;w) (2) 



Eliminating p, we find 



V-(j(r;w)~E(r;fa;) 



so that the current density can be written as 



in terms of a divergenceless vector D 



IU1 

j(r;w) = -— (D(r;u;)-E(r; W )) (3) 



V-D(r;w) = (4) 



For simplicity we assume that there is no background dielectric (es = 1) so that the polarization current is the 
mobile-charge current j. Then D has a meaning of electric displacement field for a "dielectric" medium in which all 
the current is assigned to the polarization P, namely, 

j (r; u>) = — iujP (r; u>) (5) 

so that 

D (r; w) = E (r; w) + 4ttP (r; w) (6) 

and 

V • P (r; u) = -p (r; u) (7) 

The boundary conditions reflect the facts (i) that there is diffusion so that no infinitely thin surface charge density is 
present and (ii) that there is no current through the sample boundary 

j„(r;^)| a = (8) 

Consequently, 

D n (r;w)\ g = E n (r;u)\ B (9) 

In the absence of surface charge, the electric field is continuous at the boundary 



E(r; W )| 9 = E " t (r;a;)| 9 (10) 

where E out is the field outside the particle (the sum of the applied field and the field of the particle). From eqs. ([H]) 
and (B) it then follows that the normal component of D is also continuous at the boundary. On the other hand, the 
tangential component of D may be discontinuous due to tangential currents at the boundary. 

A constitutive equation is needed to complete the description of the medium. For this we use the "generalized 
Einstein transport equation" 

j{r;cj)=a D E(r;uj)-D D Vp(r- 1 uj) (11) 

where ao is the Drude conductivity 

a D = -2_ (12) 

1 — ILOT 

ao and r are the Boltzmann conductivity and the scattering time respectively and 

Dd = —^— (13) 

1 — IWT 

where D is the static diffusion coefficient. Thus, the Einstein relation is taken to hold at w^O [E0[, namely 

a R = ao =e2 dn 

D D D ' dn K ' 

where dn/dfi is the thermodynamic density of states. 

III. SOLUTION FOR THE FIELD DISTRIBUTION 
A. General Method 

Combining eqs. (jl]), (|3|), and (|ll] ) leads to the following equation for the electric field: 

V 2 E(r; U )-A 2 E(r; U )^D(r; W ) (15) 

where A is the dynamical screening length given by 

A 2 = -^£ = A 2 (l-^) (16) 

D D \ 4tt<t d J V ' 

and 

£d = H (17) 

to 

A 2 = W^ (18) 

d[i 

are the Drude dielectric function for a conductor and the static screening length respectively. 

We introduce an auxiliary field E' 1 (r; ui) which measures the deviation from the Rayleigh-Drude field relation, 
namely, 

E /l (r;co)=E(r;^)- D(r;Cj) (19) 

(D 

It obeys the homogeneous PDE 

V 2 E h {r;cj)~A 2 E h (r-uj) =0 (20) 

and the boundary condition following from eqs. (0) and (|l0|): 



E^\ 9 =[^—^\D n \ 9 (21) 

Eha = Er t \a 

Equations (|2(j) and ( |2l|) determine the field of the screening charge completely from the boundary values of the 
displacement field D. 

Potentials can be introduced since the quasistatic fields are longitudinal: 

E = - V4>, D = - V(f> d , E h = - V(f> h (22a) 

Then from (|l|) 

cf> h = 4>-— (23) 

Given that the constant in <f> is determined, say, at infinity the undefined constants in <j> d and <j> h must compensate. 
We will see below that there is a natural choice for <f) h . The potentials are determined by solving the Laplace equation 
for (j) d in the interior: 

V 2 </> d = (24) 

V n 4> d \o = -K ut \a (25) 

plus the auxiliary PDE for <j> h : 

\7 2 cf> h -A 2 <f) h = (26) 



V„^| a = [ ^— I V„0 d | a (27) 

1 + —) \o = <f> out \o 

4> h has the special feature of determining the charge density p directly. First, from eqs. p), (0) and ( p3[ ) it follows 
that 

V 2 c/) h = -Airp (28) 

and hence from (pfl) that 

A 2 



An 



(29) 



Note that the right hand term of eq. ( |26|) is an arbitrary constant, not necessarily zero. The choice of zero, together 
with the vanishing of the total internal charge, fixes the constant in <p h to be such that the average value of <fi h over 
the sample is zero. 

B. Examples: The field in a sphere and slab 

1. Sphere 

This case is presented first because it better represents the method of potentials just described. Following Ref. Q| 
we start with the required form for the external potential of a dipole 

<r"(r;a,) = -E .r+^ (30) 

and V (|| Eo) is the induced electric dipole moment of the sphere (of radius a). The interior solution of the Laplace 
equation (E4J) for <j) d is proportional to Ep-r ex cos9 and, from the boundary condition (Ea), we get 



d (r;a;) 



-E -r 1 



2« 



(31) 



where a = P/Eq is the (complex) polarizability. Likewise <p h oc cos8. The solution of (Eu) is 

</> h (r; w) oc ii I rA j cos£? 

, . coshz sinhz 

»1 (2) = 5— 



the latter being the spherical Besscl function of the imaginary argument. The two boundary conditions (27) determine 
the coefficient of the solution <p h together with the value of a, namely, 



4> h {v-u) 



En-r 



CD-I 



3aii 



(rA) 



and 



with 



a = a 



e D + 2X sinh ( aA 



, (ep-l)-y 

6D + 2A 



coth 



X = 1-3- 



(aX) 



aA 



«A 



(32) 



(33) 



(34) 



X increases monotonically to unity as aA increases. This is very similar to the classical polarizability of a dielectric 
sphere with the dielectric constant e^ and reduces to it when the dynamical screening length is small relative to the 
radius, i.e. aA is large. The electric dipole absorption is (see Appendix A for a detailed discussion) is given by the 
imaginary part of the polarizability: 



= —uEq Im a 



The familiar RD absorption, namely, 



2„3 



Qrd = 



3u^a 

87TCTQ 



-El 



is obtained at low frequencies lu <C ctq- Corrections in (aA) are easily obtained from eq. 
order 



{RD 



11 \ 
2aA ) 



(35) 

(36) 
4|) and in the lowest 

(37) 



2. Slab 



Both the slab normal and the field are taken along z: L z <C L x ,L y ~ y/S, and —L z /2 < z < L z /2. The one- 
dimensionality permits working directly with the fields. By eqs. (Eh, (0) and ( |l0| ) D (z) = const = E out = E a 
and 



d_ 2 
dz 



-^E h (z;u))-A 2 E h (z;u>) = 



(38) 



E h (-L z /2) = E h (L z /2) 



V e D 



E 



(39) 



i c 1 \ cosh I zA ] 
EHz;u,)=E ^-V \J (40) 

£ d cosh if- A 



And the internal electric field is 



p / cosh (z A ) 



E (ziw) = -± 1 + (e D - 1) \ > (41) 

£d \ cosh%A 



The polarizability and absorption given, respectively, by 



4tt e D I L z A/2 ' 

and 

Q = 0,. (l - ^) (43) 

In summary, incorporating diffusion into the current field relation modifies the classical treatment of dielectric 
bodies by replacing the surface charge density relation a = D n (e — 1) /Aire by an additional degree of freedom for the 
bulk charge density p (r) . In both cases the Maxwell equations boundary conditions are all that is needed to assure 
self-consistency of the charge and field distribution. It is seen that the screening charge decreases exponentially with 

distance d from the surface as exp I —Ad) resulting in the correction to the RD absorption proportional to the fraction 

of the total volume occupied by the screening charge. This is a small correction for metals where A -1 is small but 
may be significant in semiconductors. 

IV. RELATION TO RESPONSE-FUNCTION FORMALISM 

In this Section we make the connection between the classical dielectric formalism in Sec. II and the quantum 
mechanical theories mentioned in the Introduction. As stated earlier, the dielectric method incorporates nonlocality 
in both the conduction and screening consistently. 

A. Response formalism 

We start with the nonlocal response function for the electric current vs. the electric field in the form first given in 
Ref. tel which was derived from a field theoretical formalism for diffusive motion, namely, 



j a (r; w)= (J aP (r, r'; u) E (r'; u) dr' (44) 

where, for both r and r' within the sample, 

a aP (r, r' ; u) = <J D (5 a p5(r - r') - V a V f3 d (r, r'; to)) (45) 

and d is the diffusion propagator (" diffuson" ) satisfying the equation 

V 2 d (r, r'; w) = -S(r - r') - ^-d (r, r'; u) (46) 

and the boundary condition 

V„d (r, r'; u) \g = V n d (r, r'; u) \ 9 > = (47) 



Outside the sample d vanishes. 

We note here two features of eqs. (44)-(fJ7|): Firstly, the conductivity tensor a with components a a p is longitudinal, 
producing longitudinal current j when E is longitudinal (= — V<^). (Furthermore, ) n = 0). Secondly, the coefficient 
with Djj on the right in eq. ( |4q ) is consistent with the appearance of ctd in eq. ( J45j ) and both are necessary. With 
them, eqs. (Ha) and (46) are the real-space equivalent of the longitudinal conductivity, eq. (3.145) of Ref. pj, in an 
unbound homogeneous medium, namely, 



o-l (q,w) 



me 



2 2 

i s q 



00 



1 — iujT + 



lD q 2 



where oq = ne 2 r /m and s is the fermi-liquid sound velocity 

2 n d/j, ne 2 Dq 

m dn m ao 



a (q, u) = a D (l - q 2 d (q, u)) 
-l 



In g-space, eqs. (W5J), (|4q) arc 



(48) 



(49) 



(50) 



d (q, u>) 



D D 



and combine to yield eq. ME 

To connect with Sec. II we next show that the constitutive relatio n (|44f) is equivalent to the Einstein transport eq. 
(PI). We need some variants of eq. (H). Inserting E = — V(j) in eq. J44|), and integrating by parts gives 



j (r; w) = iwe' 



r/// 



d(r,r / ;w)^(r / ;w)dr' 



Combining the continuity equation (g) with eq. ( |46| ) gives the nonlocal relation for the charge density 

p(r;w)= [lL{r,r';u)ct>{r';u)dr' 
where 



II (r, r'; w) = e ; 



,dn 



-5(r-r') 






d(r,r';ui) 



= e 2 ^V 2 d(r,r';a;) 



Equations ( J52J ) and (p3) give the alternative form 

Dnp (r; u) = —od4> (r; w) — iuje 2 — / d (r, r'; w) (r'; w) dr' 

the gradient of which, combined with (J5l|), yields the differential (quasilocal) constitutive ansatz (O) in Sec. II. 
The field equation ( |15[ ) is similarly verified by adding Poisson's equation to eq. (|55|), namely, 



(51) 
(52) 

(53) 
(54) 

(55) 



(V 2 -A 2 )^(r;^) = 



D D 



d(r,r';w)^(r';a;)dr' 



(56) 



This integro-differential form is typically how self-consistency appears in the microscopic theories mentioned earlier. 
But taking the gradient of (q6() (and using eqs. (J51J) and (||)) again gives a quasilocal equation which rearranges to 
(fl5|). Thus the dielectric method of Sec. II and eqs. (|44|) etc. are equivalent and the latter can be compared with the 
previous work. 

In Ref. B BM derive the unscreened polarizability of a small particle which, in the present notation, is given by 



OiOa/3 



e dn 
~2djj, 



(X a — X a ) \Xp — Xp) 



1 

V 



D D 



d(r,r';u,)s(|) 



dvdr' 



(57) 



Here, subscript "0" means unscreened and S (tH is the level density correlation function with u> — e — e' and A is 
the average level spacing. In the semiclassical limit u ~^> A, S = 1. (BM took ao to be isotropic, so their eq. (18) has 
(1/3) (r — r') ). d is d less the zero-mode term ex V~ x (see eq. (^8|) so that J ddr = 0. 

The unscreened polarizability here result from replacing <j> (r) in eq. (p2) by (f>o (r) = — r ■ Eo, with x measured from 
the "center of mass", so J x a dr — 0. The induced electric dipole is 



V Qa - / x a p (r;u>)dr = <- / XaXpTl (r, r'; u>) dv' \ E t 

— OtQafiEop 



(58) 



which is easily shown to agree with eq. (p7j) in the semiclassical limit. 

BM extended their work via the RPA towards a self-consistent theory; but their calculation were complicated and 
limited to lower order in the bubble diagrams, with the result not being self-consistent. 

A second comparison is with the papers by Mehlig and Wilkinson [g| on absorption in small particles, particularly 
the second paper in which they derive their phenomenological formalism using eqs. (E3)-(E1), but with static values 
no and Dq. Their principal result is an equation of mixed form (their eq. (9) with sign change), namely 



DV 



ioS) p = —aV 



■>AIW 



(59) 



Here 4>mw is an effective potential, taken to be the sum of the diffusive part 4> s taUc phis an ohmic part Adynamic with 
separate and pre designed frequency dependencies. Eq. ( p9[ ) was intended as the basis for a self-consistent theory. 
An equation like eq. (p9|) derives from the operation of d~ x , or V 2 + ioj/ D, on eq. (pq), giving 



V 2 - A 2 W 2 ' 



= 0= V 2 -A 2 )p 



(60) 

The right hand of (|60| ) is the time-dependent Fermi-Thomas- Debye screening equation. If A 2 is separated in its parts 
(eq. (|l6|)) and p is replaced by — V 2 (f) / Air in the first part, the form ( J59J ) results - with the differences mentioned 
above. 

Boundary conditions on cj) must be added for a complete specification of the solution of the 4th order PDE on the 
left or, correspondingly, boundary conditions on p on the right. The form of the latter is already provided in Sec. Ill 
by the boundary conditions ( |27J ) for <f> h . 

The low- frequency absorption cross-section can be obtained by expanding eq. ( p2| ) in powers of us. The result is 
that, to O (oj 2 ) , the absorption is given in terms of static quantities only and no integral equation is needed. An 
identical expression, for a special case of sphere, was given by Maksimenko et. al. in Ref. [0. Our derivation shows 
that the region of validity is uj < u>t- The details are in Appendix A. 



B. The nonlocal dielectric function 



The solution method in Sec. Ill used the auxiliary potential <fi h to bypass the need for the non-local dielectric 
function. Nonetheless, the properties of the latter are of some interest since e gives an alternative approach to 
calculating the fields and currents. The relevant dielectric tensor e (r,r';ui) is derived from the definitions dq), ra); 
and the constitutive equation is provided by eq. (|44|), giving 



D a (r; lu) = e afj (r, r'; u>) E (j (r'; u>) dv' 



with 



Airi 



e a p (r, r'; oj) = S aP S(r - r') H a af3 (r, r'; u) 



47TI. 



= e D (w) S al3 S(r - r') WaVad (r, r'; w) 



(61) 

(62) 
(63) 



The properties of a discussed above eqs. 
boundary condition (||)). 



imply that e is a longitudinal tensor (which also guaranties the 



The eigenfunction expansion of e is a natural basis for calculations. It is derived from the corresponding eigen- 
functions of d, namely, those of the Laplacian which obey the Neumann boundary conditions: eq. ( p2[ ) can be written 
as 

V 2 fl m (r) + \ 2 m n m (r) - (64) 

V n O m (r) | a = (65) 

These form a complete orthogonal set , with real and positive A 2 ^, except for the m = 0, "zero mode" 

r» = V- 1 ' 2 , A = (66) 

The eigenfunction representation of d is 

diry-.)^^ ^^ (67) 

m= m D D 

= 3(r,r'; W )+yg- (68) 

where d is the sum without m = 0. It satisfies 

V 2 d (r, r'; u) = -S(r - r ') + ±-^-d (r, v';u>) (69) 

V L>d 

and is the Fourier transform of the time-dependent diffuson of Ref . M . 
The orthonormal vector eigenfunctions are 

v m (r) = A- x VQ m (r) ; m jL (70) 

/ v m (r) • v„ (r) dr = <5„ m ; n-v m | 9 = (71) 

are complete as discussed below for the expansion of longitudinal tensor operators. 

Let A (r) = A^ (r) + At (r), longitudinal and transverse parts. The eigenfunction sum 

oo 

{ Vr,r')=E^WvM (72) 

is the longitudinal (5-function in the sense that 

6^(r,r')A (3 (r')=A La (r') 
The action of 6^g (r, r') on a longitudinal electric field is equivalent to 5 a p5 (r — r'). This leads to the representation 

oo 

e af 3 (r, r'; oj) = }^ e, n (lu) v ma (r) v mfi (r') (73) 

where 

*m M = D f--^ 6D = 1 + A 2 -^ (74) 



) D *™. - >^ Ai - -T& 



The inverse operator e is similarly expanded, with coefficients 

(e" 1 ) =e m 1 = l-A 2 l -^ (75) 

Eq. ( |75| ) may be useful for calculating the electrical field within ellipsoidal samples Q. For these cases the interior 
D-field is constant, D a —a a pEop (sum on (3), where a a g are depolarization factors. For a slab, a Q/ 3 = S a p and for a 
sphere a a/ 3 = 5 a p (l + 2V/E a 3 y Then the interior electric field is 



E a (r;u>)=^2A m (F, ;u))v ma (r) (76) 

m 

where 

A m = rrJ m • D = — J m - a -E (77) 

Cm (U) £m (w) 



and 



J m = fv m (r)dr (78) 



The frequency dependence resides entirely in the coefficients e m (w). Substituting v mQ by fi m gives the potential 
distribution 

0(r;w)=^^ m (E o ;w)n m (r) (79) 

771 

Eqs. (|7g) and ( |T9|) are not simple to evaluate except for a slab where e m and v m have simple form. The calculation 
is given in Appendix B. The solution ( ]4l| ) is reproduced. We have not carried out the sum analytically for a sphere 
and ellipsoid. 

V. ESTIMATES: DEGENERATE-ELECTRON METAL VS. NONDEGENERATE SEMICONDUCTOR 

A. Screening lengths 

Static screening length is A -1 , where for metals and semiconductors respectively, 

2 67m met e 2 
A met = — , Thomas-Fermi 

A sc = ~i — 7fT — > Debye (eo - dielectric constant) 
k B l e 



and 



Kc ("Timet 3 fc B T 

-eo 



Kaet \ n *c 2 E F 

with Timet ~ 10 22 — 10 23 cm -3 , n sc ~ 10 18 , eo ~ 10 screening length in semiconductor is ~ (10 — 100) x screening 
length in metal. For a metal, 

A mlt ~ O.lnm 

is much less than mean free path £ ~ lOnm typical for dirty metal particles. For a semiconductor, 

A" 1 - (3 - 10) nra 

It is evident that the correction 5 = 5.5/Aa in eq. (|37| ) will be too small to observe in metal particles. For a 50nm 
particle 

Smet ~ lO" 2 

On the other hand, for a semiconductor particle of the same radius 

Omet ~ .3 — 1 



10 



B. Magnetic dipole absorption 

The dipolar absorption of light is the sum of electric and magnetic parts 

1 2 1 2 

Qe = yUE Im as and Qh = x^i? Im an 

In the electric case the field is screened by the surface layer (neglecting diffusion) by a factor ~ ui/&o, while, in 
the magnetic case the currents penetrate to the skin effect depth ^/2ttoquJ/c, which is greater than the particle size. 
Taking the expressions from Ref. M §§ 72-73, we get for a sphere the ratio 

Qh_ = 8^ 2 a 2 al H 2 
Q E ~ 45 uo 2 E 2 



In free space the last factor is one. For the metallic particle (a ~ 50nm) 

Qe 

but for a semiconductor 



Qh_ 
Qe 



which can be much smaller than one. 






x above 



VI. CONCLUSIONS 



The classical Maxwell dielectric theory based on the "generalized Einstein" constitutive equation (nil) reproduces 
the results of the previous semiclassical treatments of the response of a small "dirty" conductor to an oscillating 
electric field in the limit of the continuous energy spectrum. When the discreteness of energy levels can no longer 
be neglected, the quantum theory becomes important. Both classical and semiclassical treatments require for their 
validity that the mean free path £ be small in comparison with the particle dimension, in the latter case because the 
matrix elements between impurity wave functions are treated as the excitation of the diffusion propagators. In Sec. 
IV the diffusion model, exclusive of level statistics, was shown to be equivalent to the transport equation (O). 

The real efficacy of the dielectric method is, of course, for the quantitative treatment of screening which, in both 
the nonlocal classical and quantum formalisms, involves having to solve integral equations in 3D geometries. The 
nonlocality is taken into account by the auxiliary field E^ which represents the distributions of charge (see eq. (p9[)). 
Solving for the fields is not significantly more complicated than for bodies with local dielectric constants. 

An unanswered question in the above considerations is the validity of the screening description at short distances. 
There are three lengths of interest: the mean free path £, the screening distance A -1 and the particle size a. In 
metallic particles (£ <IC a) typically A -1 is much less than t and the transport equation (O) has questionable physical 
meaning on the scale of A -1 . Nonetheless, in the limit A -1 — > where all charge is on the surface, the classical 
Rayleigh-Drude theory of metallic particle is obtained. This suggests that the correction for finite screening in eqs. 
(]43|) and (J37 ) is still qualitatively correct. The same criticism applies to the quantum mechanical descriptions in Refs. 
Pp H and [8j as well as to Pines-Nozieres equation (Eq). The theory here may provide a more consistent picture for 
semiconductor particles where screening lengths are much larger. 

It is interesting that, despite the above mentioned physical limitation, the mean free path itself does not appear 
explicitly — only scattering time r in the Drude form 1 — ilot. Nor does the Thouless frequency ujt appear explicitly. 
It does determine the applicability of the low-frequency absorption in Appendix A. 

Finally we conjecture that in the extreme quantum region when the levels are discrete and lu and/or T < A, use 
of the "Golden Rule" and impurity-averaged spectral quantities like S (w/A) is fundamentally wrong. We will return 
to this point in later work. 



If 



APPENDIX A: LOW-FREQUENCY ABSORPTION 

First, we review expressions for absorption of energy from an applied electric field. The local absorption rate is 

g(r; W ) = iR e {j(r; W ).E*(r; W )} (Al) 

The integrated absorption rate 

Q(w)= J q(r;uj)dr (A2) 

can be rewritten for a longitudinal electric field by applying Gauss' theorem and the continuity equation, 

Q(cj)=^lmU p(r;u,)4>*(r;uj)dr\ (A3) 



The symmetry of this form is used in the derivation of eq. (All) below. 

A less symmetric form for Q (us) in which cf> (r) is replaced by c/)q (r), the applied field in the absence of the sample, 
is familiar from thermodynamic arguments (Ref. M) for quasistatic fields. We remark that replacing <j) by cf>o is valid 
at all frequencies because Coulomb forces are conservative and do no net work in a cycle. We digress briefly to prove 
this (probably well-known) statement, namely that the induced field 

4>ind = 0-00 

does no work. Because the system is isolated electrically, the total induced charge pi n d is zero and 4>i n d ~ r~ 2 at large 
distances. (An initial charge distribution po docs not contribute to Q and is ignored.) So pi n( i = p and 

V 2 0= VVmd = spinet (A4) 

with different conditions at "infinity." To see that 

AQ = / Im {pindfand} = (A5) 

integrate by parts over all space. The surface term goes to zero at large r, and the new integrand Im < |V0j n( j| \ 

vanishes. 

For the small systems considered here <pQ (r) = — Eo-r so that 

QH = |lmp.E*} (A6) 

where V is the induced moment V = J rpi„d (r; u>) dr, which is J P (r; u>) dr for the dielectric medium (eq. (ffl)). For 
linear response, P=aEo, ( |A6| ) becomes 

QH = ||E | 2 Im{aM} (A7) 

To lowest order in u, Q oc oj 2 . A symmetrical form of the coefficient of uj 2 can be derived from eq. ( |A3| ) by 
expanding p and 4> to first order in to. Write 

p = p s t+upi, (j> = (j) s t +u><f>i (A8) 

Now <j) st = <f) h because cd (0) = oo, so that 

A 2 

Pat = —r<t>«t ( A9 ) 

47T 



The product of the static terms is real so it gives zero in (|A3|), as expected, and leaves 

Q^ImjyV^-p^drj 



(A10) 
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To 0(u) eqs. (52)-(64[) give 



where 



and 



Pi = J [n st (r, r') cf >1 (r') + U x (r, r') <f> at (r')] dv' 



U st (r, r') = n (r, r'; 0) = -e 2 ^S (r - r') 



(All) 



(A12) 



U 1 (r,r')=- — — d(r,r';0) 
D d\i 



(A13) 



The term linear in <fii in the brackets of eq. (All) is real so that the absorption, to order w 2 , is determined by u) = 
quantities 



ui 



A 



<P* st (r)d(r,r';0)<p st (r')dr'dr 



where 



e" an \ f 2 dn\ 1 A 4 



D dfi (Tq V dn 



Co (47r) 



(A14) 



(A15) 



From the eigenfunction expansion of d(r,r';ui), eq. (pd), we see that d(u>) sw d(0) as long as lo < DX\ ~ Wy, 
the Thouless frequency. The Drude frequency correction wt to ctq & n d £* is negligible at these frequencies since 



lu t t w {(■/ L) . The l atter must be small in order for the constitutive equation (|ll|) to be meaningful physically. 



We have checked (A14) for the spherical case (eq. (|33|)). From eq. (p2|) 



rt (r) =(-Eo-p)/(r) 

/ (r) = — csch (aA) i\ (rA) 



(A16) 
(A17) 



d (r, r'; 0) is the seen from eqs. (pig), ([47]) be the Coulomb Green function with Neumann boundary conditions 



d(r,r';0) = - 



1 



4-7T r — r 



1 /r 2 

2V \!$ 



E 



« + 1) (rrj 

A 1 2/ + 1 P ' C0S< 



where P; (cos0) is the Lagrange polynomial and 9 is the angle between r and r'. Only the I = 1 component 



di(r,r';0) 



l/lr< 



D V47TT?, 



27ra c 



cos 6* 



(A18) 



(A19) 



(where r< and r> denote the lesser and greater, respectively, of r and r') contributes to the integral of eq. ( A14| ). 
The result 



,2„3 



Q = |^|E | 2 F(A a ) 

87TCT0 

= Q RD F (Aa) 



where 



F (z) = 1 ( coth z ) + ^r ( coth z ) H coth z (coth 2 z - l) 

z \ z J z z \ z J 2z 2 ' 



(A20) 

(A21) 

(A22) 



agrees with the term Im j_o| (B3I) which is linear in ui. 

Finally, we compare (A14) with the expression for absorption obtained by Maksimenko et. al.. in Rcf. [Q who 
extended the semiclassical method of GE [pi to take screening into account. The unscreened dipolar energy — er • Eq 
is replaced by a "screened" dipole — eR • Eq. The latter is found within the RPA in which the irreducible polarization 
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part has the same form as polarization itself, namely, having matrix elements between disordered states. The solution 
of the RPA integral equatio n, af te r an a verage over states at the Fermi surface and restricting to spherical geometry, 



has the same form as eqs. (A16), (A17): 



<t> st (r) =-R(r)-E (A23) 

R(r)=r/(r) (A24) 

The agreement reflects the correspondence between the self-consistency of the RPA and of the Maxwell dielectric 
theory. (The resulting formula for Q in Ref. [H contains the level correlation function S (w/A) and reduces to the 
present form when S — 1.) 

APPENDIX B: DIELECTRIC TENSOR IN A SLAB 

The eigenfunctions and eigenvalues of d(r,r';u) are Corresponding to eqs. ( p4[) and 

«(«,/;«) = y^cos \—lL\--< z <- (Bl) 



and 



The vector eigenfunctions are 



■ 2 . W / T17T \ 2 . U) 



X --^ = {-) ^DZ ^ 



v m (z)=z^-sin(^— U+2JJ ( B3 ) 

We treat the case in Sec. Ill where Eo||z so the vector notation z can be dropped. Then D = Eo and the prescription 
in eq. (|7|) gives 

fe^ = 1 V -J_ ( 1 - ^ 1 , ^ _1 sin f ^±11 ) ( B4) 

^ n^2n+l\ ((2n+l)7r/L) 2 +A 2 ; \ I / 

The summation is easily done and yields the solution (M 
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